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1. Introduction 

De Broglie wave interferometry using cold atoms is emerging as a new tool for basic 
science and technology. There are numerous approaches and applications which have 
evolved since the first demonstration experiments in the early 1990's. These notes will 
not attempt an exhaustive or comprehensive survey of the field. An excellent overview is 
provided in Ref. [1] and other lectures in this volume. These notes will focus on what has 
become known as light-pulse atom interferometry, which has found fruitful applications 
in gravitational physics and inertial sensor development. 

These notes are organized as follows. We we first summarize basic theoretical con- 
cepts. We will then illustrate this formalism with a discussion of applications in inertial 
navigation and in a detailed design discussion of an experiment currently under develop- 
ment to test the Weak Equivalence Principle. 

2. Atom interferometry overview 

By analogy with their optical counterparts, atom interferometery seeks to exploit 
wave interference as a tool for precision metrology. In this measurement paradigm, 
a single particle (photon or atom) is made to coherently propagate over two paths. 
These paths are subsequently recombined using a beamsplitter, and their relative phase 
becomes manifest in the probability of detecting the particle in a given output port of 
the device. Hence, measuring the particle flux at the interferometer output ports enables 
determination of the phase shift. As this relative phase depends on physical interactions 
over the propagation paths, this enables the characterization of these interactions. 
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A key challenge for de Broglie wave interferometers is the development of techniques 
to coherently divide and recombine atomic wavepackcts. For simplicity, consider an 
atom with initial momentum p, characterized by a wavefunction tp(x) ~ exp \ip ■ x] 
The atom then is subject to a Hamiltonian interaction which is engineered to evolve 
the wavepacket into a momentum superposition state. One such interaction spatially 
modulates the amplitude of the wavefunction, so that, for example, exp [j-p • x] — ► 
/(x) exp \ip ■ xl , where /(x) is a real periodic function with spatial frequency k. Fourier 
decomposing /(x) immediately shows that the final wavefunction is a coherent super- 
position of momenta p, p ± ftk, p ± 2ftk 7 etc. In practice, such an interaction can be 
implemented by passing a collimated atomic beam through a microfabricated transmis- 
sion grating, as demonstrated by Pritchard and co-workers [2] . 

Another interaction is one which spatially modulates the phase of the wavefunction. 
Consider, for example, ex P [*k ' X ]V>( X )- This interaction results in a momentum 

translation p — > p + Kk. A particularly useful implementation of this process imparts 
a spatial phase modulation by driving transitions between internal atomic states. For 
simplicity, consider a two-level atom with internal states |1) and |2) that are resonantly 
coupled by an applied optical traveling wave E oc exp [ik ■ x] via the clcctric-dipole 
interaction /t-E (where /t is the dipolc moment operator). If the atom is initially prepared 
in state ?A( X ) |1)> then following an interaction time t its state becomes a(f)?/>(x) |1) + 
(3{t) exp [ik • x]^>(x) |2) (see Section 3 for details). The interaction time can be chosen, 
for example, so that \a\ — = 1/V2 to implement a beamsplitter (the tt/2 pulse 
condition). In this case, the internal state of the atom becomes correlated with its 
external momentum. In practice, two-photon stimulated Raman transitions between 
groundstate hypcrfinc levels have proven to be particularly fruitful for implementing this 
class of beamsplitter. Why? Transitions are made between long lived hypcrfine levels 
while the phase grating periodicity is twice that of a single photon optical transition 
(when the Raman transition is driven in a counter-propagating beam geometry) . 

The above mechanisms operate in free space. A new family of atom optics, based on 
control of atom wavepacket motion in atomic waveguides, is under development. The 
basic idea is that atoms are steered using microfabricated wires deposited on surfaces. 
These are loosely analogous to optical fiber waveguides for light. In principle, coherent 
beamsplitters are implemented by the appropriate joining of waveguides. Recently, a 
combination of microwave and magnetic fields has allowed for creation of waveguide 
structures capable of coherent wavefront division. These structures have been used to 
demonstrate proof-of-principle interferometer topologies which have been used to study 
the coherence properties of the division process. The notes below will discuss free space 
atom optics, which have proven effective for precision measurement of incrtial forces. 

Exploiting the momentum exchange principles outlined above, it becomes straight- 
forward to devise a de Broglie wave interferometer which is based on sequences of light 



( x ) In reality, the wavefunction is spatially localized, so by the uncertainty principle, there must 
be a corresponding spread in momentum about the mean value. 



4 



Jason M. Hogan, David M. S. Johnson and Mark A. Kasevich 



pulses. For example, consider a three pulse sequence based on Raman transitions. An 
initial Raman tt/2 pulse places an atom in a coherent superposition of wavepackets in 
states |1) and |2) whose mean momenta differ by ?ik ff (here k e g is the effective wavevec- 
tor of the Raman process, see below). After an interrogation time T these wavepackets 
separate by a distance hk c sT/m, where m is the atomic mass. A subsequent optical pulse 
is then applied whose duration is chosen to drive the transitions |1) — > |2) and \2) — > |1) 
with unit probability (a it pulse). This pulse has the effect of redirecting the momenta 
of the wavepackets so that at a time T later the wavepackets again overlap. A final n/2 
pulse then serves as the exit beamsplitter. 

3. Phase shift determination 

In this section we review the method for calculating the phase difference between the 
two halves of the atom at the end of the light-pulse atom interferometer pulse sequence 
outlined above. These results are well-known [3, 4], but we are not aware of a complete, 
formal derivation of these rules in the literature. Other equivalent formalisms for this 
calculation do exist (see, for example [5, 6]). For Section 5'2 it is necessary to understand 
the formulae for the phase difference (Section 3T). The proof of these formulae as well 
as a discussion of their range of validity is given in Section 3'2 but is not necessary for 
the rest of the paper. 

3T. Phase shift formulae. - The main result we will show is that the total phase 
difference A(j) to t between the two paths of an atom interferometer may be written as the 
sum of three easily calculated components: 

(1) A</>tot = A</> propaga t 

ion ~t~ A<^> SC p ara tj on -\- A(^>i aser . 

For this calculation we take h = c = 1. 

The propagation phase A<^ propaga tion arises from the free-fall evolution of the atom 
between light pulses and is given by 

(2) A^ propagation = ( [ ( L c- E ^ dt ) ~ 2 ( / {Lc - E ^ dt 

upper ^"'fr ' lower ^ tj 

where the sums are over all the path segments of the upper and lower arms of the 
interferometer, and L c is the classical Lagrangian evaluated along the classical trajectory 
of each path segment. In addition to the classical action, Eq. (2) includes a contribution 
from the internal atomic energy level £7j. The initial and final times tj and tF for each 
path segment, as well as L c and E i7 all depend on the path segment. 

The laser phase A</>i ascr comes from the interaction of the atom with the laser field 
used to manipulate the wavefunction at each of the beamsplitters and mirrors in the 
interferometer. At each interaction point, the component of the state that changes 
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momentum due to the light acquires the phase of the laser </>l(£o> x c(*o)) = k • x c (t ) 
ujt n + <j> evaluated at the classical point of the interaction: 



(3) A0 lascr = \J2± ( f> L (t j ,xL u (t j )) - \J2±ct> L (t j , Xl (t j )) 



lower 



The sums are over all the interaction points at the times tj , and x u (t) and x; (t) are the 
classical trajectories of the upper and lower arm of the interferometer, respectively. The 
sign of each term depends on whether the atom gains (+) or loses (— ) momentum as a 
result of the interaction. 

The separation phase A0 scpara ti O n arises when the classical trajectories of the two 
arms of the interferometer do not exactly intersect at the final beamsplitter (see Fig. 1). 
For a separation between the upper and lower arms of Ax = x; — x„, the resulting phase 
shift is 

(4) A</> SO p aration = p • Ax 

where p is the average classical canonical momentum of the atom after the final beam- 
splitter. 

3'2. Justification of phase shift formulae. - The interferometer calculation amounts 
to solving the Schrodinger equation with the following Hamiltonian: 

(5) #tot = H a + H cxt + V int (x) 

Here H a is the internal atomic structure Hamiltonian, H cx t is the Hamiltonian for 
the atom's external degrees of freedom (center of mass position and momentum), and 
Vint(x) = — /t • E(x) is the atom-light interaction, which we take to be the electric dipole 
Hamiltonian with fj, the dipole moment operator. 

The calculation is naturally divided into a series of light pulses during which V\ n t ^ 0, 
and the segments between light pulses during which Vint = and the atom is in free-fall. 
When the light is off, the atom's internal and external degrees of freedom are decoupled. 
The internal eigenstates satisfy 

(6) id t \Ai) = H a \Ai) = Ei \Ai) 

and we write the solution as \Ai) = \i) e - lB «(*-*o) w itfi time-independent eigenstate |z) 
and energy level 

For the external state solution we initially consider H cxt — H(x,p) to be an 
arbitrary function of the external position and momentum operators: 



(7) 



idt\1>)=H(i,p)\il>) 



6 



Jason M. Hogan, David M. S. Johnson and Mark A. Kasevich 



It is now useful to introduce a Galilean transformation operator 

(8) G c = G(x c , Pc ,L c ) = e i/ic* e -*-xc e ip e -* 

which consists of momentum boost by p c , a position translation by x c , and a phase shift. 
We choose to write 

(9) |V> = Gc \<Pcm) ■ 

We will show that for a large class of relevant Hamiltonians, if x c , p c , and L c are taken 
to be the classical position, momentum and Lagrangian, respectively, then \4>cm) is a 
wavepacket with (x) = (p) =0, and the dynamics of \4>cm) do not affect the phase 
shift result (i.e., \4>cm) is the center of mass frame wavefunction) . However, for now 
we maintain generality and just treat x c , p c , and L c as arbitrary functions of time. 
Combining (7) and (9) results in 

(10) id t \<Pcm) = {g\H(x, p)G c - i& c dtG c } \4>cm) 

= { H ( x + X c , p + p c ) + p c • X - (p + p c ) • X c + L c } \(j>CM) 

where we used the following identities: 

(11) G c xG c = x + x c 

G c pG c = p + p c 

G\H{x, p)G c = H(x + x c , p + Pc) 

Next, we Taylor expand H(x + x c , p + p c ) about x c and p c , 

(12) H(x + x c , p + p c ) = H(x c , p c ) + (x c , p c ) • x + Vpi?(x c , Pc) • P + H 2 

where 77 2 contains all terms that are second order or higher in x and p. (We will 
ultimately be allowed to neglect H2 in this calculation.) Inserting this expansion and 
grouping terms yields 

id t \4>cm) = I (H c - x c • p c + L c ) + (V Xc H c + p c ) • x + (V Pc II c - x c ) • p + H 2 } \4>cm) 

where we have defined the classical Hamiltonian H c = H(x c , p c ). If we now let x c , p c , 
and L c satisfy Hamilton's equations, 



(13) x c 

Pc 



= V pc #c 
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with p c = V± a L c the classical canonical momentum, then \4>cm) must satisfy 
(14) id t \4> C m) = H 2 \4> CM ) 

Next we show that it is possible to choose \4>cm) with (x) = (p) = for a certain 
class of H 2 , so that x c and p c completely describe the atom's classical center of mass 
trajectory. This is known as the semi-classical limit. Starting from Ehrcnfcst's theorem 
for the expectation values of \4>cm), 



(15) d t (xi) = 1 

(16) d t (pi) = 1 

and expanding about (x) and (p), 



H 2 ,Xi 
H 2 ,pi 



'pi 

= - (d Xi H 2 



d t (x^ = ( dp i H 2 



<*>,<p> 



+ d Pj d Pt H 2 



(Pj - (Pj)) + d X] d Pt H 2 



dt (pi) 



+ 



d Xi H 2 
1 



2! 



<*>,<p> 
d Xk d Xj d Xi H 2 



(i),(p> 

d Xi d Xi H 2 



<*>,<P> 

(Pj ~ <Pj» (Pk ~ <Pfc» + • • 



<*>,<p> 



(% - (%)) 



(xj - (xj)) + dp d Xz H 2 



<i).<P> 



<*},<P> 

% - (%>) (£fc - + 



<x),<P> 



we find the following: 
(17) d t (xt) = d p M 2 

W,(p) 

(is) a, = - d Xi H 2 



21 



dp k dp j dp z H 2 



<*>,<p) 



... ... ~ — ■ <>■> Hi 

WAP) ^! 



<*>,<P> 



Ax% + ■■■ 



where Ax^ k = (xjXk) — (%) (£fc) and Ap 2 - k = (pjpk) — (Pj) (Pk) & re measures of the 
wavepacket's width in phase space ( 2 ). This shows that if H 2 contains no terms higher 
than second order in x and p, then Ehrcnfcst's theorem reduces to Hamilton's equations, 
and the expectation values follow the classical trajectories. Furthermore, this implies that 
we can choose \</>cm) to be the wavefunction in the atom's rest frame, since (x) = (p) = 
is a valid solution to Eqs. (17) and (18) so long as all derivatives of H 2 higher than 
second order vanish. In addition, even when this condition is not strictly met, it is 
often possible to ignore the non-classical corrections to the trajectory so long as the 
phase space widths Axjk and Apjk are small compared to the relevant derivatives of 
H 2 (i.e., the semi-classical approximation). For example, such corrections are present 



( ) In general, there will also be cross terms with phase space width such as (xjpk) — {pk)- 
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for an atom propagating in the non-uniform gravitational field g of the Earth for which 
dfdfdfH2 ~ d 2 g. Assuming an atom wavepacket width Ax < 1 mm, the deviation 
from the classical trajectory is d t {p} ~ (d^g)Ax 2 < 10~ 20 g, which is a negligibly small 
correction even in the context of the ~ 10~ 15 g apparatus we describe below for testing 
the Equivalence Principle. 

The complete solution for the external wavefunction requires a solution of Eq. (14) 
for \4>cm), but this is non-trivial for general Hi- In the simplified case where Hi is 
second order in x and p, the exact expression for the propagator is known [7] and may 
be used to determine the phase acquired by \4>cm)- However, this step is not necessary 
for our purpose, because for second order external Hamiltonians the operator H2 does 
not depend on either x c or p c . In this restricted case, the solution for the rest frame 
wavefunction \4>cm) does not depend on the atom's trajectory. Therefore, any additional 
phase evolution in \4>cm) must be the same for both arms of the interferometer and so 
does not contribute to the phase difference. This argument breaks down for more general 
H2, as does the semi-classical description of the atom's motion, but the corrections will 
depend on the width of \4>cm) m phase space as shown in Eqs. (17) and (18). We ignore 
all such wavepacket-structure induced phase shifts in this analysis by assuming that the 
relevant moments {Axjk, Apjk, ■ ■ .} are sufficiently small so that these corrections can 
be neglected. As shown above for the non-uniform (d 2 g ^ 0) gravitational field of the 
Earth, this condition is easily met in many experimentally relevant situations. 

Finally, we can write the complete solution for the free propagation between the light 
pulses: 

(19) (x\il>,Ai) = (x|G c |0 CM > |^) = e^^yP-^-^cMtx-Xc) \i) e -™^-u) 

We see that this result takes the form of a traveling wave with de Broglie wavelength 
set by p c multiplied by an envelope function <J>cm (x), both of which move along the 
classical path x c . Also, the wavepacket accumulates a propagation phase shift given by 
the classical action along this path, as well as an additional phase shift arising from the 
internal atomic energy: 

(20) A^ propagation = ( [ ( L c- E ^ dt ) - S ( / {Lc ~ E ^ dt ) 

upper ^"'fr ' lower ^ tj ' 

where the sums are over all the path segments of the upper and lower arms of the 
interferometer, and tj, tF, L c , and Ei all depend on the path. 

Next, we consider the time evolution while the light is on and Vi n t 7^ 0. In this case, 
the atom's internal and external degrees of freedom are coupled by the electric dipolc 
interaction, so we work in the interaction picture using the following state ansatz: 



(21) 



|*) = jdvY,Ci(\>)\i>v)\Ai) 
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where we have used the momentum space representation of \()>cm) and so \ip p ) = G c e lH2 ^ 
Inserting this state into the Schrodinger equation gives the interaction picture equations, 

(22) id t |*) = tjrfp^ d -^- |V P ) \A) + H & |*) + H cxt |*> = H tot |*> 

i 

(23) => Ci(p) ee ^ = 7 / dp ' E C ^P') (^l (V'pI ^t(x) | V> |^-> 



where we used (6) and (7) as well as the orthonormality of \Ai) and \ip P }- The interaction 
matrix element can be further simplified by substituting in \ip p ) and using identity (11): 

(24) (Vp|Mnt(x) IV) = (p|e i ^(*- t °)M nt (x + x c ) e -^(*-*o) |p') 

, if J^-S^H (t-t ) 

= <p|v5 nt (x+x c )| P ') e V 2 ™ 2 ™; 1 

where we have made the simplifying approximation that Hi J^. This approximation 
works well as long as the light pulse time t = t — t is short compared to the time 
scale associated with the terms dropped from Hi- For example, for an atom in the 
gravitational field of Earth, this approximation ignores the contribution m(d r g)x 2 from 
the gravity gradient, which for an atom of size Ax 1 mm leads to a frequency shift 
~ m(d r g)Ax 2 ~ 1 mHz. For a typical pulse time t < 1 ms, the resulting errors are 
< 1 /xrad and can usually be neglected. Generally, in this analysis we will assume the 
short pulse (small r) limit and ignore all effects that depend on the finite length of the 
light pulse. These systematic effects can sometimes be important, but they are calculated 
elsewhere [8] [9] . In the case of the 87 Rb- 85 Rb Equivalence Principle experiment we discuss 
below, such errors are common-mode suppressed in the differential signal because we use 
the same laser pulse to manipulate both atoms simultaneously. 

As mentioned before, we typically use a two photon process for the atom optics 
(i.e., Raman or Bragg) in order to avoid transferring population to the short-lived ex- 
cited state. However, from the point of view of the current analysis, these three-level 
systems can typically be reduced to effective two-level systems [10] [11]. Since the result- 
ing phase shift rules are identical, we will assume a two-level atom coupled to a single 
laser frequency to simplify the analysis. Assuming a single traveling wave excitation 
E(x) = E cos (k • x - ut + <j>), Eq. (23) becomes 
(25) 

*(P) = | /rfp'E ^ c ^P') <Pl (e^'^)-^) + h.c) |p') e l ^ + ^'^ dt 

3 

where the Rabi frequency is defined as fiy = (i\ (— /t • Eq) \j) and u>ij = Ei — Ej. Now 
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we insert the identity 

(26) k • (x + x c ) - tot + <f> = k • x + (k ■ x c (t ) - ut Q + + (k • x c - u)dt 

boost laser phase Doppler shift 

into Eq. (25) and perform the integration over p' using (p| e ±lk x |p') = (p|p' ± k): 

(27) c,(p) = 1^ SUj { C ,(p-k) e ^e l ^-" +k -^ + ^-^ )dt + 

■jip + k) e - i ^ e -^^-^- w+k -^ + ^+^) d * J 



j 

< 



where we define the laser phase at point {t ,x c (to)} as (j>L = k • x c (t ) — ut^ + (f>. Finally, 
we impose the two-level constraint (i = 1,2) and consider the coupling between ci(p) 
and c 2 (p + k) : 

(28) ci(p) = -nc 2 (p + k)e- i * i e~ <J * A(p)<it 

2i 

c 2 (p + k) = ■ifi*c 1 (p)e i *V / * A(p)<tt 
2i 

Here the detuning is A(p) = lo — to + k • (x c + ^) + the Rabi frequency is Q = 
Q12 = (f2 2 i)*, and io$ = LO21 > 0. In arriving at Eqs. (28) we made the rotating wave 
approximation [12], dropping terms that oscillate at (u>o+u>) compared to those oscillating 
at (lo — to). Also, flu = since the \Ai) are eigenstates of parity and /t is odd. 

The general solution to (28) is 
(29) 

ci(p,i) \ = ( A c (p) e -^ A (P)- - t A s (p)e-^ A ^e-^\ / Cl ( Pjto ) 
v c 2 (p + k,f)y { -iA:(p)e* A (p) T e^ A*(p)e2 A (p)- J ^c 2 (p + k,t ) 



(30) A c (p) = cos( ^ A(p)i + t ) sinl -VA(p) 2 + |^| 2 r 



\^ { py + M \) +i ^m 

> V A (P) 2 + 1 01 " 

(31) A s (p) = " sin (iv/A(p) 2 + |^| 2 r 

^/A(p) 2 + |Q| 2 ^ 

In integrating (28) we applied the short pulse limit in the sense of k • x c r 2 <C 1, ignoring 
changes of the atom's velocity during the pulse. For an atom falling in the gravitational 
field of the Earth, even for pulse times r ~ 10 ^s this term is ~ kgr 2 <~ 10~ 2 rad which 
is non-negligible at our level of required precision. However, for pedagogical reasons we 
ignore this error here. Corrections due to the finite pulse time are suppressed in the 
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proposed differential measurement between Rb isotopes since we use the same laser to 
simultaneously manipulate both species (see Section 5'1). 

For simplicity, from now on we assume the light pulses are on resonance: A(0) = 0. 
We also take the short pulse limit in the sense of |A(p) — A(0)| r <C 1 so that we can 
ignore all detuning systematics. This condition is automatically satisfied experimentally, 
since only the momentum states that fall within the Dopplcr width <~ t _1 of the pulse 
will interact efficiently with the light. 
(32) 

ci(p,t) \ _f Ac -iA s e-^\ ( ci(p,to) "\ A c = cos^ 



c 2 (p + k,t); \-iA*e^ A c J \c 2 (p + k,t )) A, = -fasm^ 

In the case of a beamsplitter (J pulse), we choose \Q\t = §, whereas for a mirror (it 
pulse) we set \Q \ t = n: 

(33) A/ -f * f^ L \ A -( ° ~ ie " 

These matrices encode the rules for the imprinting of the laser's phase on the atom: the 
component of the atom Ci(p, to) that gains momentum from the light (absorbs a photon) 
picks up a phase +<Pl, and the component of the atom C2(p + k, to) that loses momentum 
to the light (emits a photon) picks up a phase — <f>^. Symbolically, 

(34) |p)^|p + k)e^ 

(35) |p + k) -» \p)e~^ L 

As a result, the total laser phase shift is 

(36) A0 lascr = ( J2±4> L (t j ,x u (t j )) ] - ( ^±0 L (t,,x ; (t,)) 

lower 

where the sums are overall of the atom-laser interaction points {tj,x u (tj)} and {tj,xi(tj)} 
along the upper and lower arms, respectively, and the sign is determined by Eqs. (34) 
and (35). 

The final contribution to A0 to t is the separation phase, A</> sep aration- As shown in 
Fig. 1 , this shift arises because the endpoints of the two arms of the interferometer need 
not coincide at the time of the final beamsplitter. To derive the expression for separation 
phase, we write the state of the atom at time t = t + t just after the final beamsplitter 
pulse as 



(37) 



|*(t)> = |*„(t)> + !*,(*)> 
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Fig. 1. - Separation Phase. This is a magnified view of the end of the interferometer which 
shows the upper and lower arms converging at the final beamsplitter at time to, and the result- 
ing interference. The dashed and solid lines designate the components of the wavefunction in 
internal states \Ai) and 1^-2), respectively. After the beamsplitter, each output port consists of 
a superposition of wavepackets from the upper and lower arm. Any offset Ax = x; — x u between 
the centers of the wavepacket contributions to a given output port results in a separation phase 
shift. 



where \^ u (t)) and \^i(t)} are the components of the final state that originate from the 
upper and lower arms, respectively. Just before the final beamsplitter pulse is applied, 
we write the state of each arm as 

(38) |*»(*o)> = / <feci(p,to)Gu(to) |p) |^i)e ie " 

(39) |*j(*o))= /"dpc 2 (p,fo)G(to)|p)|^2)e ie! 

where G u = G(x„, p u , L u ) and G/ = G(x/, pi, Li) are the Galilean transformation opera- 
tors for the upper and lower arm, respectively. These operators translate each wavepacket 
in phase space to the appropriate position (x u or xj) and momentum (p u or p;). Here 
we have assumed for clarity that prior to the final beamsplitter the upper and lower arms 
are in internal states \A\) and \A%) with amplitudes ci(p, to) and C2(p, to), respectively; 
identical results are obtained in the reversed case. We have also explicitly factored out 
the dynamical phases 9 U and 9i accumulated along the upper and lower arms, respec- 
tively, which contain by definition all contributions to laser phase and propagation phase 
acquired prior to the final beamsplitter. 
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We write the wavefunction components after the beamsplitter in the form of Eq. (21): 

(40) |tt u (i)) = [dpJ2 #\p,t)Gu |P) \Ai) 

(41) |*,(*)> = /dpE c i°(P'*)^IP>l^> 

where we invoked the short pulse limit so that e~ lH2T w 1. Next we time evolve the 
states using Eq. (32) assuming a perfect -| pulse and using the initial conditions given in 
Eqs. (38) and (39): namely, c^\p,t a ) — ci(p,t o )e I0u and c^\p, to) = for the upper 
arm and c^(p,t ) = and c%\p,to) — C2(p,t )e ldl for the lower arm. 

(42) |* u (t)> = |rfp C i(p,i )|^G„|p>|yl 1 ) + ^ e ^( x ")G u |p + k)|A 2 )|e^ 

(43) |tt,(t)> = |dp C2 (p + k,i ){^e-^( x ')G ; |p)|^i) + ^G ; |P + k)|A 2 )|e^ 

We now project into position space and perform the p integrals, 



(44) (x|*„(t)) = Cl(x Jt uM) { e <p «-( x - x ») - i e *M*«) e i(p»+k)-(*-*») 

(45) (x|*,(i)) = C2(X ~J^ M) |- ie -^( x ')e l (P'- k )'( x - x ') |Ax) + e*P<'( x - x ') |A 2 )} e l,?! 

where we identified the Fourier transformed amplitudes using Cj(x— x c , to) = /^P ( x — x c|p) Cj(p, to). 
The resulting interference pattern in position space is therefore 

(x|*(t)) = <x|* u (t)) + <x|* I (t)) 

= \A Y ) {c 1 (x-x„,t )e l6 '"e 4p "-( x - x ") - l c 2 (x-x 1 ,t )e ai e-^ (x ' ) e , ( [ "- k H x - x ')} 

+ -J= |A 2 ) {c 2 (x - xi, io)e^e lp! - (x - Xi) - i ci(x - x„, f ) e «- e i * i(x - ) e i(p " +k) - (x - x - ) } 
The probability of finding the atom in either output port \A\) or |A 2 ) is 

(46) |<^| (x|*(i)>| 2 = |Cl|2 + |C2 ' 2 + \ (id + c.c.) 

(47) |<A 2 | (x|*(t))| 2 = |Cl|2 + |C2 ' 2 - 1 (» Cl ^e^ +c.c.) 
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with Ci = ci(x — x„, to) and C2 = c 2 (x — x;, i )- For the total phase shift we find 
(48) A0! = \0 U + p„ • (x - x u )} - {0, - </> L (x,) + (p, - k) • (x - x ; )} 



(49) = 8 u -[6i- fo(xz) + pi • Ax + Ap • (x - x 



A<^p r0 p aga ti on , 1 + A0iascr,l A0 SO paration, 1 

and 

(50) A0 2 = {6 u + <j)L{*u) + (Pu+k) • (x-x„)} - \0i + Pl • (x-x;)} 

(51) = (0 U + <f> L (xL u )^ -61 + p 2 -Ax + Ap-(x-x) 

^.^propagation, 2 + A0 lasC r,2 A0 sop a r atio„,2 

where pi = p "+( P '~ k ) and P2 = ( p " +k ) +p ' arc the average momenta in the \A\) (slow) 
and 1^2) (fast) output ports, respectively. In general, the separation phase is 

(52) A<^ S eparation = P • Ax 

which depends on the separation Ax = x; — x u between the centers of the wavepackets 
from each arm as well as the average canonical momentum p in the output port. We 
point out that even though the definitions (48) and (50) use the same sign convention 
as our previous expressions for laser (36) and propagation (20) phase in the sense of 
( )u PP or - ( ) lower, the separation vector Ax is defined as (x) lowcr - (x) uppcr . 

Notice that the phase shift expressions (49) and (51) contain a position dependent 
piece Ap • (x - x), where x = x "+ x ' and Ap = (p„ + k) - p ; = p u - (p; - k), 
owing to the fact that the contributions from each arm may have different momenta 
after the last beamsplitter. Typically this momentum difference is very small, so the 
resulting phase variation has a wavelength that is large compared to the spatial extent 
of the wavefunction. Furthermore, this effect vanishes completely in the case of spatially 
averaged detection over a symmetric wavefunction. 

Finally, we show that the total phase shifts A^i and Atfo for the two output ports 
are actually equal, as required by conservation of probability. According to Eqs. (49) 
and (51), the contributions to the total phase differ in the following ways: 

^A(/>p ropa g a tion,l + A01 aseri i^ — ^A0p ropa g a tion,2 + A(/> iasor ,2^ = <Al,( x /) — </>l(x m ) 

= k • (xj - x„) = k • Ax 

Acceptation, 1 - A(/> scparat ion,2 = Pi ' Ax - p 2 ■ Ax = -k • Ax 

Together these results imply that A<pi = Afe and prove that the total interferometer 
phase shift At/> t ot is independent of the output port. 
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The accuracy of the above formalism is dependent on the applicability of the aforemen- 
tioned stationary phase approximation as well as the short pulse limit. The stationary 
phase approximation breaks down when the external Hamiltonian varies rapidly com- 
pared to the phase space width of the atom wavepacket. The short pulse limit requires 
that the atom's velocity not change appreciably during the duration of the atom-light 
interaction. Both approximations are justified to a large degree for a typical light pulse 
atom interferometer, but in the most extreme high precision applications such as we con- 
sider here, important corrections are present. However, we emphasize that these errors 
due to finite pulse duration and wavepacket size are well-known, previously established 
backgrounds. 

4. Applications in inertial navigation 

The navigation problem is easily stated: How do we determine a platform's trajec- 
tory as a function of time? In the 20th century, solutions to the problem have led to the 
development of exquisitely refined hardware systems and navigation algorithms. Today 
we take for granted that a hand-held GPS reciever can be used to obtain meter level 
position determination. When GPS is unavailable (for example, when satellites are not 
in view), position determination becomes much less accurate. In this case, stand-alone 
"black-box" inertial navigation systems, comprised of a suite of gyroscope and acclerom- 
eters, are used to infer position changes by integrating the outputs of these inertial force 
sensors. State-of-the-art commercial grade navigation systems have position drift errors 
of kilometers per hour of navigation time, significantly worse than the GPS solution. 
How can we close the performance gap between GPS and inertial systems? The way 
forward is improved instrumentation: better gyroscopes and accelerometers. 

The light-pulse interferometry method is well suited to inertial applications. As shown 
in Section 3T, the phase shift for the light pulse interferometer consists of contributions 
from path phases, optical interactions, and separation phases. However, for the sensitivity 
range of interest, the optical phase shifts dominate, and therefore A0 tot w A(j>i aser . In 
this case there is a straightforward interpretation for the operation of the sensor: the 
sensor registers the time evolution of the relative position of the atomic wavepackets with 
respect to the sensor case (defined by the opto-mechanical hardware for the laser beams) 
using optical interferometry. Since distances are measured in terms of the wavelength of 
light, and the atom is in a benign environment (spurious forces, such as from magnetic 
field gradients, are extremely small - below I0~ w g), the sensors are characterized by 
superbly stable, low-noise operation. 

The phase shift between the two paths is inferred by measuring the probability of 
finding the atoms in a given output port. Since we are only concerned with the optical 
phase shift, we can calculate the sensor output using Eq. (3). We see that for a standard 
7r/2 — 7r — 7r/2 excitation sequence [35], A0; aser = ki • (x" — xj) — k 2 • (x£ — x°) — 
k 2 • (x 2 — x 2 ) + k 3 • (x 3 — x°) . Here the subscript indexes each of the three successive 
optical interactions, x" and x' are the semi-classical positions of the atom along the 
upper and lower paths, respectively, at the time of each interaction in a non-rotating, 
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Fig. 2. - Schematic illustration of an atomic beam gyroscope. 



inertial coordinate system, x° is the phase reference for the optical fields, and is the 
propagation vector of the laser field associated with each pulse. 

4T. Gyroscope. - Assuming the atoms have initial velocity vo and the effective Raman 
propagation vectors initially have common orientations k which rotate with angular rate 
^2, it is straightforward to show that A(j)i aser w k • (2v x fl)T 2 . This configuration 
is well-suited to precision measurements of platform rotations. This expression can be 
put in a form analogous to the Sagnac shift for optical interferometers by noting that for 
small rotation rates 0, A = ^(k x v )T 2 is the area enclosed by the interfering paths. 
Thus this shift can also be written in the Sagnac form <f> = ^pfi • A - proportional to 
the product of the enclosed area and rotation rate. 

Gyroscopes built on this principle have achieved performance levels in the laboratory 
which compare favorably with state-of-the art gyroscopes, as shown in Figs. 2, 3, 4 and 5. 
Key figures of merit for gyroscope performance include gyroscope noise (often referred to 
as angle random walk), bias stability (stability of output for a null input) and scale factor 
stability (stability of the multiplier between the input rotation rate and output phase 
shift). The laboratory gyroscope illustrated below has achieved a demonstrated angle 
random walk of 3 ^(deg/hr 1 ^ 2 , bias stability of 60 ^(deg/hr (upper limit) and scale factor 
stability of 5 ppm (upper limit). Key drivers in the stability of the gyroscope outputs 
are the stability of the intensities used to drive the Raman transitions, and alignment 
stability of the Raman beam optical paths. Non-inertial phase shifts associated with 
magnetic field inhomogeneities and spurious AC Stark shifts are nulled using a case 
reversal technical where the propagation directions of the Raman beams are periodically 
reversed using electo-optic methods. 

4'2. Accelerometer . - If the platform containing the laser beams accelerates, or if 
the atoms are subject to a gravitational acceleration, the laser phase shift then contains 
acceleration terms A<p = k • aT 2 . For a stationary interferometer, with the laser beams 
vertically directed, this phase shift measures the acceleration due to gravity. Remarkably, 
part per billion level agreement has been demonstrated between the output of an atom 
interferometer gravimeter and a conventional, "falling-corner-cube" gravimeter. [14]. In 
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the future, compact, geophysical (10~ 8 g accuracy) grade instruments should enable low 
cost gravity field surveys (see Fig. 6). For this type of instrument, laser cooling meth- 
ods are used to initially prepare ensembles of roughly 10 7 atoms at kinetic temperatures 
approaching 1 /iK. At these low temperatures, the rms velocity spread of the atomic 
ensembles is a few cm/sec. Cold atom ensembles are then launched on ballistic trajecto- 
ries. In this configuration, the time between laser pulses exceeds 100 msec, which means 
wavepackets separate by roughly 1 mm over the course of the interferometer sequence. 
The phase shift is read-out by detecting the number of atoms in each final state using 
resonance fluorescence and normalized detection methods. 

In general, both rotation and acceleration terms are present in the sensor outputs. 
For navigation applications, the rotation response needs to be isolated from the accel- 
eration response. In practice, this is accomplished by using multiple atom sources and 
laser beam propagation axes. For example, for the gyroscope illustrated in Fig. 2, 
counter-propagating atom beams are used to isolate rotation induced phase shifts from 
acceleration induced shifts. It is interesting to note that the same apparatus is capable 
of simultaneous rotation and acceleration outputs - a significant benefit for navigation 
applications which require simultaneous output of rotation rate and acceleration for three 
mutually orthogonal axes. Since gyroscope and accelerometer operation rest on common 
principles and common hardware implementations, integration of sensors into a full in- 
ertial base is straightforward. Of course, particular hardware implementations depend 
on the navigation platform (e.g. ship, plane, land vehicle) and trajectory dynamics. 
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Fig. 4. - Gyroscope interference fringes (a) and power spectral density of gyroscope noise (b). 
The increase in noise in the 10-50 Hz band is due to technical rotation noise sources in the 
building (the gyroscope was mounted directly to the building floor). See Ref. [13]. 



4'3. Gravity gradiometer. - There is an additional complication in navigation system 
architecture for high accuracy navigation applications: the so-called "problem of the 
vertical." Terrestrial navigation requires determining platform position in the gravity 
field of the Earth. Due to the Equivalence Principle, navigation system accelerometers 
do not distinguish between the acceleration due to gravity and platform acceleration. So 
in order to determine platform trajectory in an Earth-fixed coordinate system, the local 
acceleration due to gravity needs to be subtracted from accelerometer output in order to 
determine the acceleration of the vehicle with respect to the Earth. This means that the 
local acceleration due to gravity needs to be independently known. For example, existing 
navigation systems use a gravity map to make this compensation. However, in present 
systems, this map does not have enough resolution or accuracy for meter-level position 
determination. To give a feeling for orders of magnitude, a 1CT 7 error in knowledge of 
the local acceleration due to gravity integrates to meter-level position errors in 1 hour. 

There are at least two paths forward: 1) better maps or 2) on-the-fly gravity field 
determination. Improved maps can be obtained with more precise surveys. On-the- 
fly determination seems impossible, due to the Equivalence Principle (since platform 
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Fig. 5. - (a) Drift in case reversed gyroscope output with time and its correlation with instrument 
temperature (solid black), (b) Gyro noise output residual after subtracting temperature model, 
(c) Allan deviation of gyro output before and after temperature compensation (see Ref. [13]). 



accelerations cannot be discriminated from the acceleration due to gravity). However, 
the outputs from a gravity gradiometer - an instrument which measures changes in the 
acceleration due to gravity over fixed baselines - can be used for this purpose. The 
idea is to integrate the gravity gradient over the inferred trajectory to determine gravity 
as a function of position. In principle, such an instrument can function on a moving 
platform, since platform accelerations cancel as a common mode when the output from 
spatially separated accelerometers are differenced to obtain the gradient. In practice, such 
a strategy places hard requirements on the stability of the component accelerometers: 
their responses need to be matched to an exceptional degree in order to discriminate 
gravity gradient induced accelerations (typically below 10~ 9 g) from other sensor error 
sources. 

Due to the stability of their acceleration outputs, a pair of light pulse accelerometers 
is well-suited to gravity gradient instrumentation. The basic idea is to simultaneously 
create two spatially separated interferometers using a common laser beam. In this way, 
technical acceleration noise of the measurement platform is a common-mode noise source 
which leads to near identical phase shifts in each accelerometer. On the other hand, a 
gravity gradient along the measurement axis results in a residual differential phase shift. 
This configuration has been used to measure the gravity gradient of the Earth, as well 
as the gravity gradient associated with nearby mass distributions, as illustrated in Fig. 
7. Laboratory gravity gradiometers have achieved resolutions below 1 E (where 1 E = 
1CP 9 scc~ 2 ). This configuration has also been used to measure the Newtonian constant 
of gravity G [15, 16, 17]. In future navigation systems, an ensemble of accelerometers, 
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configured along independent measurement axes, could acquire the full gravity gradient 
tensor. 

In addition to their role as navigation aids, gravity gradiometers have applications 
in geodesy and oil/mineral exploration. The idea here is that mass/density anomalies 
associated with interesting geophysical features (such as kimberlite pipes - in the case of 
diamond exploration - or salt domes - in the case of oil exploration) manifest as gravity 
anomalies. In some cases, these anomalies can be pronounced enough to be detected 
by a gravity gradiometer from an airborne platform. Atom-based gravity gradiometers 
appear to have competitive performance figures of merit for these applications as com- 
pared with existing technologies. For these applications, the central design challenge is 
realization of an instrument which has very good noise performance, but also is capable 
of sustained operation on a moving platform. Figs. 8, 9 illustrate a system currently 
under development for this purpose. 
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Fig. 7. - Vertical axis gravity gradiometer. Two atom interferometer accelerometers separated 
along a vertical measurement axis are simultaneously interrogated by a common set of Raman 
laser beams. This apparatus was used to measure G, Newton's constant [15]. 




Fig. 8. - Gravity gradiometer with horizontal measurement axis. Raman beams propagate along 
the axis defined by the tube connecting each accelerometer housing. 
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Fig. 9. - Gravity gradiometer response to a proof mass which is periodically brought in close 
proximity to one of the accelerometer regions. For the proof masses used for this demonstration, 
the apparatus is capable of resolving changes in G at the 3 x 1CT 4 level [18]. 



5. Application to tests of the Equivalence Principle 

Precision tests of the Equivalence Principle (EP) promise to provide insight into 
fundamental physics. Since the EP is one of the central axioms of general relativity 
(GR), these experiments are powerful checks of gravity and can tightly constrain new 
theories. Furthermore, EP experiments test for hypothetical fifth forces since many 
examples of new forces are EP-violating[19]. 

The Equivalence Principle has several forms, with varying degrees of universality. 
Here we consider tests of the Weak Equivalence Principle, which can be stated as fol- 
lows: the motion of a body in a gravitational field in any local region of space-time is 
indistinguishable from its motion in a uniformly accelerated frame. This implies that the 
body's inertial mass is equal to its gravitational mass, and that all bodies locally fall at 
the same rate under gravity, independent of their mass or composition. 

The results of EP experiments are typically expressed in terms of the Eotvos param- 
eter r\ = Aa/a, where Aa is the EP violating differential acceleration between the two 
test bodies and a is their average acceleration [20]. Currently, two conceptually differ- 
ent experiments set the best limits on the Equivalence Principle. Lunar Laser Ranging 
(LLR), which tests the EP by comparing the acceleration of the Earth and Moon as they 
fall toward the Sun, limits EP violation at 7] = (-1.0 ± 1.4) x 10~ 13 [21]. Recently, the 
Eot-Wash group has set a limit of rj = (0.3 ± 1.8) x 10~ 13 using an Earth-based torsion 
pendulum apparatus [22] . Several proposed satellite missions aim to improve on these lim- 
its by observing the motion of macroscopic test bodies in orbit around the Earth [23, 24]. 
Here we discuss our effort to perform a ground-based EP test using individual atoms 
with a goal of measuring r\ <~ 10~ 15 . Instead of macroscopic test masses, we compare 
the simultaneous acceleration under gravity of freely-falling cold atom clouds of 87 Rb and 
85 Rb using light-pulse atom interferometry [25]. 

Light-pulse atom interferometers have already been used to make extremely accu- 
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rate inertial force measurements in a variety of configurations, including gyroscopes, 
gradiometers, and gravimeters. For example, the local gravitational acceleration g of 
freely- falling Cs atoms was measured with an accuracy ~ 10~ 9 g [14]. Gravity gradiome- 
ters have been used to suppress noise as well as many systematic errors that are present 
in absolute g measurements by comparing the acceleration of two displaced samples of 
atoms. A differential measurement of this kind was used to measure the Newtonian 
constant of gravity G with an accuracy of <~ 3 x 10~ 3 G [15]. The EP measurement 
we describe here benefits from an analogous differential measurements strategy, where 
in this case the common-mode noise suppression arises from a comparison between two 
co-located isotopes of different mass, rather than between spatially separated atoms as 
in a traditional gradiometer. 

5T. Proposed experiment overview. - The proposed experiment ideally consists of 
simultaneously observing the free-fall motion of the two Rb isotopes in the absence of 
all non-gravitational forces. To this end, the measurement is performed inside a 10 cm 
diameter by 8.8 m long cylindrical ultra high vacuum chamber. To maximize their free- 
fall time, the atoms are launched in a vertical fountain geometry from the bottom of 
the chamber. Light-pulse atom interferometry is performed while the atoms are in free- 
fall, and the resulting phase shift is sensitive to the atoms' acceleration. Figure 10 is a 
space-time diagram depicting the trajectories that each atom follows during the free-fall 
interferometry sequence. 

To maximize the cancellation of spurious effects, both the 87 Rb and 85 Rb atom clouds 
are launched at the same time and are made to follow the same trajectories as closely as 
possible. We launch both isotopes from the same magnetic trap in order to minimize any 
differences between their initial positions and velocities. As a result of the small isotope 
shift between 87 Rb and 85 Rb, we are able to use the same laser pulses to simultaneously 
manipulate them during the interferometer sequence (see Fig. 10). Using the same laser 
makes the apparatus insensitive to pulse timing jitter and dramatically reduces the phase 
noise stability requirements of the lasers. 

A single measurement of acceleration in our atom interferometer consists of three 
steps: atom cloud preparation, interferometer pulse sequence, and detection. In the 
first step a sub-microkelvin cloud of ~ 10 7 atoms is formed using laser cooling and 
evaporative cooling in a TOP trap [26]. This dilute ensemble of cold atoms is then 
launched vertically with velocity v z <~ 13 m/s by transferring momentum from laser 
light using an accelerated optical lattice potential [27]. This technique allows for precise 
control of the launch velocity and, because it is a coherent process, it avoids heating the 
cloud via spontaneous emission. 

In the second phase of the measurement, the atoms follow free-fall trajectories and 
the interferometry is performed. A sequence of laser pulses serve as beamsplitters and 
mirrors that coherently divide each atom's wavepacket and then later recombine it to 
produce interference (see Fig. 10). The atom beamsplitter is typically implemented us- 
ing a stimulated two-photon process (Raman or Bragg transitions), resulting in a net 
momentum transfer of k e ff = k2 — ki w 2k2 at each interaction. Since the acceleration 
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Fig. 10. - Space-time diagram for a single atom of mass m during the interferometer pulse 
sequence. The atom is launched with velocity v z from the bottom of the vacuum system. At 
time t = 0, a -| (beamsplitter) pulse is applied to coherently divide the atom wavefunction. 
After a time T, a n (mirror) pulse is applied that reverses the relative velocity between the 
wavefunction components. A final ^ (beamsplitter) pulse at time 2T results in interference 
between the two space-time paths. The interferometer phase shift is inferred by measuring the 
probability of detecting the atom in either state |1) (solid line) or state 2) (dashed line). Note 
that points D and E are in general spatially separated in the presence of non-uniform forces, 
leading to a separation phase shift. 



sensitivity of the interferometer is proportional to the effective momentum hk e g trans- 
ferred to the atom during interactions with the laser, we intend to take advantage of 
more sophisticated atom optics. Large momentum transfer (LMT) beamsplitters with 
hk c ff — 24hk have been demonstrated [28], and up to lOOftfc may be possible. Promis- 
ing LMT beamsplitter candidates include optical lattice manipulations [27] , sequences of 
Raman pulses [29] and adiabatic passage methods [30]. 

The third and final step of each acceleration measurement is atom detection. At the 
end of the interferometer sequence, each atom is in a superposition of the two output 
velocity states, as shown by the diverging paths on the right in Fig. 10. These two final 
velocity states are directly analogous to the two output ports of a Mach-Zehnder light 
interferometer after the final recombining beamsplitter. As with a light interferometer, 
the probability that an atom will be found in a particular output port depends on the 
relative phase acquired along the two paths of the atom interferometer. Since the output 
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states differ in velocity by ~ Hk e s/m, they spatially separate over time. After an appro- 
priate drift time, the two velocity groups can be separately resolved, and the populations 
can be then measured by fluorescence imaging. 

We now consider the expected sensitivity of our differential 87 Rb- 85 Rb accelerometer. 
Recent atom interferometers have demonstrated sensor noise levels limited only by the 
quantum projection noise of the atoms (atom shot noise) [31]. Assuming a time-average 
atom flux of n — 10 6 atoms/s, the resulting shot noise-limited phase sensitivity is ~ 
= 1CP 3 rad/\/Hz. The phase shift in an atom interferometer due to a constant 
acceleration g is A<f> — k c sgT 2 [4]. Taking advantage of the L w 8.8 m vacuum system 
allows for a long interrogation time of up to T = 1.34 s. Finally, using hk e s — lOftfc 
LMT beamsplitters results in an acceleration sensitivity of <~ 7 x 10~ 13 g/VHz and 
a precision of < lQ~ 15 g after ~ 1 day of integration. In the most conservative case, 
constraining ourselves to conventional 2hk atom optics leads to a precision of < 10~ 15 g 
after ~ 1 month of integration. This estimate is based on realistic extrapolations from 
current performance levels, which are at 10~ 10 g [15]. 

5'2. Error model. - An accurate test of the EP requires a thorough understanding 
of potential backgrounds. To reach the goal sensitivity, we must control spurious ac- 
celerations to < 10~ 15 g. Systematic errors at this level can arise from many sources, 
including gravity gradients, Earth's rotation, and electromagnetic forces. To calculate 
these contributions to the phase shift, we follow the prescription outlined in Section 3. 
We take the atom's Lagrangian in the lab frame to be 

(53) L = lm(r + ft x (r + R e )) 2 - m0(r + R e ) - i«B(r) 2 

where r is the position of the atom in the lab frame, R e — (0, 0, R e ) is the radius of 
the Earth, ft = (0,Q, y ,Q z ) is the Earth's rotation rate, and cf)(r) is the gravitational 
potential. In the chosen coordinate system, z is the vertical direction in the lab and ft 
lies in the y-z plane. We then expand cf) in a Taylor scries about R e , 

(54) <f>(r + R e ) = - ^g • r + ^(Ty)^- + ^(Qijk)nrjr k + ^(%fei)^r J r fcn 

where Earth's gravity field is g = — V0(R e ), the gravity gradient tensor is = d^gu the 
second gradient tensor is Qijk = dkdjgi, the third gradient tensor is Sijki = didkdjgi, and 
repeated indices are summed over. Since z is in the vertical direction in the lab we have 
that g = (0,0, —g) and g = 9.8 m/s 2 . The interferometer follows a fountain geometry 
which is approximately one-dimensional along the z-direction, so we only include Qzzz 
and S zzzz and safely ignore the other second and third gradient tensor terms. Likewise, 
in this analysis we assume that off-diagonal gradient tensor terms with i ^ j are 
small and can be ignored (this is exactly true for a perfectly spherical Earth). The 
effects of higher-order moments of the gravitational field are treated separately using a 
perturbative calculation as described in Section 5'2.1. 
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Because magnetic fields can cause significant systematic errors, the atoms are prepared 
in one of the magnetic field insensitive clock states (\rriF = 0) states). The residual energy 
shift in a magnetic field B is then Ub = |aB 2 , where a is the second order Zeeman shift 
coefficient. We consider magnetic fields of the form 

(55) B(r) = (b + ^-z) z 



dz 

where B is a constant bias magnetic field, and ^ is the gradient of the background 
magnetic field. While this linear model is sufficient for slowly varying fields, in Sec- 
tion 5'2.2 we describe a perturbative calculation that can account for more complicated 
magnetic field spatial profiles. 

We do not include additional electromagnetic forces in the Lagrangian as their accel- 
erations are well below our systematic threshold. For neutral atoms, electric fields are 
generally not a concern since the atom's response is second order. Furthermore, electric 
fields are easily screened by the metallic vacuum chamber, leading to negligibly small 
phase shifts. Short range effects due to the Casimir [32] force or local patch potentials 
[33] are also negligible since the atoms are kept far (> 1 cm) from all surfaces throughout 
the experiment. 

As explained in Section 3'2, we point out that the phase shift derived from Eq. 1 is 
only exactly correct for Lagrangians that are second order in position and velocity. When 
this is not true, as is the case in Eq. 53 when Qijk ^ and Sijki ^ 0, the semiclassical 
formalism breaks down and there are quantum corrections to the phase shift. However, 
these corrections depend on the size of the atom wavepacket compared to the length scale 
of variation of the potential, and are typically negligible for wavepackets ~ 1 mm in size. 

To analytically determine the trajectories r(t), we solve the Euler-Lagrange equations 
using a power series expansion in t: 

N 

(56) ri(t) = J2ain(t-t ) n (i = 1,2,3) 

The coefficients a,i n are determined recursively after substitution into the equations of 



g dr n 



< 1. For 



motion. This expansion converges quickly as long as flT <C 1 and 
our apparatus with characteristic length r ~ 10 m and time T ~ 1 s these conditions 
are easily met, since QT <~ 10~ 4 rad and y-§^ir ~ (ir) ~ assuming a spherical 
Earth. With these trajectories and the interferometer geometry shown in Fig. 10 we 
obtain the following expressions for the phase shift in the slow (state |1)) output port: 

(57) A0p ropa g a tion = jX(Sac + Sce) — (Sab + Sbd)) 

(58) A</> lascr = <f> L (r A , 0) - L (r c , T) - <f> L (r B , T) + 4>l (r D , 2T) 

(59) A(/> scparation = 7^(Pz3 + Pe) ■ (r D - r E ) 
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where Sij is the classical action along the path segment between points r, and r 3 , and 
Pi = diL(ri) is the classical canonical momentum at point after the final beamsplitter. 
The laser phase shift at each interaction point is 

(60) cj) L (r, t) = k eff • r - u) ctf t + O 

where k e ff and uj c s are the effective propagation vector and frequency, respectively, for 
whatever atom-laser interaction is used to implement the atom optics. In the case of the 
stimulated two-photon processes mentioned earlier, k e ff = k 2 — ki and w e s = u>2 — u)\ = 
(k 2 - h)/c. 

Using the above method, we computed the phase shift response for a single atom 
interferometer, and the results are shown in Table I. The values of the experimental 
parameters used to generate this list are representative of the 8.8 m apparatus described 
previously. Many of these terms are common to both species, and in order obtain our 
< 10~ 15 g sensitivity, we rely on their common mode cancellation. In Table II we compute 
the differential phase shift between a 87 Rb and a 85 Rb interferometer. The two species 
have different masses m and second order Zeeman coefficients a, as well as potentially 
different launch kinematics r(0) and f(0). To create Table II, we parameterized the 
launch kinematics with a differential velocity 5v — (5v x ,Sv y ,Sv z ) and initial position 
Sr = (Sx, 5y, Sz) between the centroids of the two isotope clouds. Residual systematic 
phase errors are the result of differential accelerations that arise from gravity gradients, 
second gravity gradients, coriolis and centrifugal forces, and magnetic forces on the atoms. 

As justified below, we expect to achieve experimental parameters that reduce the 
majority of the systematic errors below our experimental threshold. However, the first 
several terms in Table II are still too large. In order to further reduce these backgrounds, 
we can employ propagation reversal to suppress all terms cx fc^ ff . This well-known tech- 
nique entails reversing the laser propagation vector k c ff — ► — k e ff on subsequent trials 
and then subtracting the two results [15]. This suppresses terms 1, 4, 9, and 10 by 
Ak e ff/k e s, where Afc e ff is the error in fc e ff made as a result of the reversal. Reducing 
these terms below our systematic threshold requires Ak e g/k e s < 10~ 5 . The main accel- 
eration signal and all other terms linear in k e g are not suppressed by this subtraction. 

After propagation vector reversal, the last important background phase shifts arise 
from the differential coriolis and centrifugal acceleration between the isotopes (Table II 
terms 2, 5, 7, and 8), and from the Earth's gravity gradient (Table II terms 3 and 6). 
We discuss the techniques used to control these remaining systematics in Section 5'3. 

5'2.1. Gravity Inhomogeneities. The Taylor series expansion of the gravitational 
potential (see Eq. 54) is a good approximation of the coarse structure of Earth's gravity 
on length scales of ~ R e , the radius of the Earth. However, local gravity can also vary 
on much shorter length scales in a way that depends on the specific mass distribution 
surrounding the experiment, and these gravity inhomogeneities can result in spurious 
phase shifts. Since these inhomogeneities can be rapidly spatially varying, the Taylor 
series expansion is not well-suited for their description. Instead, we leverage the fact 
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Phase shift Size (rad) Fractional size 



1 


-k cS gT 2 


-2.85 x 10 8 


1.00 


2 


kefiRetiyT 2 


6.18 x 10 5 


2.17 x 10" 3 


3 


-k cH T zz v z T 3 


1.58 x 10 3 


5.54 x 10 -6 


4 


J2 keagTzzT 4 


-9.21 x 10 2 


3.23 x 10 -6 


5 


-3k cfi v z Q,lT 3 


-5.14 


1.80 x 10 -8 


6 


2k e gv x Q y T 2 


3.35 


1.18 x 10 -8 


7 


lk Bsg a 2 y T 4 


3.00 


1.05 x 10 -8 


8 


1~2 ^ c ff-^ e -^~' 22: 


2.00 


7.01 x 10" 9 


9 


ft^eff T" 1 T'S 


7.05 x 10 _1 


2.48 x 10" 9 


10 


f keffgQzzzVzT 5 


9.84 x 10 -3 


3.46 x 10" 11 


11 


1^^°^ Q z z zV Z T 
4 &eff -^e ^ 


-7.66 x 10" 3 


2.69 x 10" 11 


12 


-6.50 x 10" 3 


2.28 x 10" 11 


13 


— Y20^effp QzzzT 


-3.81 x 10" 3 


1.34 x 10" 11 


14 


-3.39 x 10" 3 


1.19 x 10" 11 


15 


2m y 


-2.30 x 10 -3 


8.06 x 10" 12 


16 


^k e sT zz v z T 


2.19 x 10" 3 


7.68 x 10" 12 


17 


— 360 kcftgT zz T 


-7.53 x 10" 4 


2.65 x 10" 12 


18 


3k c gVyQ. y Q z T 3 


2.98 x 10" 4 


1.05 x 10" 12 


1 Q 


b „0 O 1/ T 2 


7 a 1 v 1 n -5 


9 fin v 1 n -13 

Z.OU A 1U 


on 


3 £• „B D ,, n 2 T 5 

4 ft'cff zzz U z i.Ly± 


u y 1 n -5 


7 Kfi v 1 (T 14 

J .JU A ±U 


9 1 


31 £- „nB D 2 T 6 


1 /i 7 v 1 rr 5 


t 17 v 1 n -14 

0.1/ A 1U 


22 


-kcftT^,vM 2 T'' 

9 "'Cli ZZ V "y 


-1.42 x 10" 5 


5.00 x 10" 14 


23 


"QfctffTzz'VxQyT 


1.08 x 10 -5 


3.81 x 10" 14 


24 




-6.92 x 10" 6 


2.43 x 10" 14 


25 


7hk cif (-\ ». 7^4 
? 12m 4 


-6.84 x 10 -6 


2.40 x 10" 14 


26 


g keffTxxVxQy T 


-5.42 x 10 -6 


1.90 x 10" 14 


27 


-§k cS gTzzn 2 yT 6 


4.90 x 10" 6 


1.72 x 10" 14 


28 




4.75 x 10" 6 


1.67 x 10" 14 


29 


3hk cfi ^ 1 ^-,5 
8m IJ^zzzJ- 


4.40 x 10" 6 


1.55 x 10" 14 


30 


fceff i?e T" 2 z i7 y T 
eff T - ^ t~>5 


1.63 x 10" 6 


5.74 x 10" 15 


31 


-1.63 x 10 -6 


5.74 x 10" 15 


32 


9.78 x 10" 7 


3.43 x 10" 15 


33 


lifc cff aB (^-B)T 2 


-7.67 x 10" 8 


2.69 x 10" 16 


34 


QQfceffQ$zzzz'V z T 


-7.52 x 10 -8 


2.64 x 10" 16 


35 


4 &eff 5j zzzV z T 


3.64 x 10" 8 


1.28 x 10" 16 


36 


j2 keff T ZZ Q zzz v Z T 


-3.13 x 10 -8 


1.10 x 10" 16 



Table I. - Phase shift response for a single atom interferometer f — tt — f sequence given 
the Lagrangian in Eq. 53. Column 3 shows the fractional size of each term compared to the 
acceleration signal k e jjgT 2 . All terms with fractional phase shift > 10~ 16 are included. The 
numbers are for a 87 Rb interferometer with the following parameters: k e g — 2k — 2 • 7S 27T lm , 
T zz = -2g/R e , T xx = T yy = g/R e , Q zzz = Gg/R 2 e , S Z zzz = ~24g/Rl, R e = 6.72 x lQ S m, 
Bo = 100 nT, and d z B = 0.1 nT/m. The Earth's rotation rate is given by £l y — flcosdLat and 
Q.z = Osin^Lat with — 7.27 x 10~ 5 rad/ s and 0Lat = 37.4 degrees North latitude. The initial 
position of the atom in the lab is taken as r(0) = (xo,yo:0), with x — 1 mm and yo — 1 mm. 
The initial velocity is f(0) = (v x ,v y ,v z ), with v x = 1 mm/ s, v y — 1 mm/ s, and v z = 13.2 m/s. 
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Phase shift Size (rad) Fractional size 



2 2k e s8v x Q y T 

3 -/c off T^5w z T 3 

4 -§f^ ^) hk 2 eS Q 2 y T 3 

2 ^ m s5 m 87 y en a 

5 -3fc e ff0 2 <5i; z T 3 

6 -k cB T zz 8zT 2 

7 -k cS 8yQ. y Q. z T 2 

8 3k eS Sv y n y Q z T 3 

11 -(it-^)^ (a z B)r 2 

12 —2k c gT xx SxQ y T 



3 



1.66 x 10" 2 


5.83 


x 


lO" 11 


3.35 x irr 3 


1.18 


x 


lO" 11 


1.44 x 10" 4 


5.05 


X 


lO" 12 


-5.40 x 10" 5 


1.90 


X 


lO" 13 


-4.68 x 10 -6 


1.64 


X 


lO" 14 


8.93 x 10" 7 


3.14 


X 


lO" 15 


-7.41 x 10" 7 


2.60 


X 


1Q -15 


2.98 x 10" 7 


1.05 


X 


10 -is 


-1.61 x 10" 7 


5.65 


X 


1Q -16 


1.03 x 10" 7 


3.63 


X 


1Q -16 


-9.94 x 10" 8 


3.49 


X 


1Q -16 


-6.92 x 10" 8 


2.43 


X 


1Q -16 



Table II. - Differential phase shift between 87 Rb and 85 Rb. To create the differential phase shift 
list we parameterized the launch kinematics with a differential velocity (8v x = 1 fim/s, 8v y = 
1 \xml s, 8v z — 12 fim/s) and position (Sx = 1 /im, 8y = 1 fim, 8z — 10 nm) between the centroids 
of the two isotope clouds. All other parameters are the same as in Table I. Column 3 shows the 
fractional size of each term compared to the acceleration signal k e ggT 2 . We include all terms 



with a fractional phase shift > 10 



that these inhomogcncitics arc typically small in magnitude and solve for the induced 
phase shift using first-order perturbation theory [3] . This linearization allows us to make 
a Fourier decomposition of the phase shift response in terms of the spatial wavelengths 
of the local g-field. 

First, we assume a one-dimensional gravitational potential perturbation of the form 
5(j)(z). The gravity field perturbation along the vertical (z) direction is defined as Sg z (z) = 
—d z S4> and may be written as 

(61) Sg z (z) - J Sg z (X)e^dX 

where Sg z (X) is the Fourier component of a gravity perturbation with wavelength A. The 
total phase shift due to gravity inhomogeneities summed over all wavelengths is 

(62) A^ g = J T gz (X)Sg z (X)d\ 

where T gz (X) is the interferometer's gravity perturbation response function. Qualita- 
tively, the response to short wavelengths is suppressed since the interferometer averages 
over variations that are smaller than its length [34] . The response is flat for wavelengths 
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Fig. 11. - Differential gravity response function versus spatial wavelength A between simulta- 
neous 87 Rb and 85 Rb interferometers. Short wavelengths are averaged over by each individual 
interferometer, while long wavelength inhomogeneities cancel as a common-mode between the 
two species. This response curve assumes identical launch kinematics for the two isotopes. 



longer than the scale of the interferometer, and in the limit where A ~ R e this analysis 
smoothly approaches the results of our Taylor series calculation described above. 

For the 87 Rb- 85 Rb EP measurement, we are interested in the differential phase re- 
sponse between the isotopes. Figure 11 shows the differential response function AT gz (X) = 
\(T gz )s7 — {T gz )85\ for gravity inhomogeneities. Once again, short wavelength variations 
are suppressed since each interferometer spatially averages over a ~ 10 m region. The 
peak response occurs at a length scale set by the spatial separation of the arms of a 
single interferometer Az = MrftX' ^ 16 mm. Perfect differential cancellation between 
isotopes is not achieved because the spatial separation of the arms is mass dependent. 
Additionally, the long wavelength differential response is suppressed because the differ- 
ences between the isotope trajectories are negligible when compared to variations with 
length scales much longer than Az. 

The differential response curve allows us to compute systematic errors arising from the 
specific gravity environment of our interferometer. Quantitative estimates of these effects 
requires knowledge of the local 5g z (z), which may be obtained through a combination of 
modelling and characterization. The atom interferometer itself can be used as a precision 
gravimeter for mapping Sg z (z) in situ. By varying the launch velocity, initial vertical 
position, and interrogation time T, the position of each gravity measurement can be 
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A [meters] 

Fig. 12. - The magnitude power spectra of the local gravitational field, 8g z , for several example 
sources. The solid (blue) curve is a 10 -2 kg point source, 10 cm from the center of the inter- 
ferometer. Similarly, the dotted (purple) curve is a 1 kg source at 1 m and the dash-dotted 
(yellow) curve is 1000 kg at 10 m. The long-dashed (green) curve is a thin 10 m long rod of 
mass 10 kg, parallel to the interferometer, whose center is 1 m from the interferometer. 



controlled. 

Figure 11 shows that the differential 87 Rb- 85 Rb interferometer is maximally sensitive 
to short wavelength (A ~ 1 — 10 cm) gravitational inhomogeneities. To investigate the 
impact of local uneven mass distributions on the experiment, we compute the spectrum 
Sg z (X) of various sources at different distances from the interferometer. These results 
are shown in Fig. 12. When combined with our response function (Fig. 11), we see that 
for typical mass inhomogeneities, only those that are within a few centimeters of the 
interferometer can cause potentially significant systematic phase shifts. These nearby 
inhomogeneities result in phase errors of ~ 10~ 6 rad, which is slightly above our target 
sensitivity. It will therefore be especially important for the EP measurement that we 
characterize the local g-field at the centimeter scale. 

5 '2. 2. Magnetic field inhomogeneities. The linear expansion of B in Eq. 55 approx- 
imates large scale variation of the magnetic field. However, local field inhomogeneities 
may exist on short length scales due to the presence of nearby magnetic materials. These 
variations are not well approximated by a Taylor series expansion. Using the same pro- 
cedure described above for gravity inhomogeneities, we write the local magnetic field 
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Fig. 13. - Differential magnetic field response function between simultaneous Rb and Rb 
interferometers. Short wavelengths are averaged over by each individual interferometer, while 
long wavelength inhomogeneities cancel as a common-mode between the two species. This 
response curve assumes identical launch kinematics for the two isotopes. 



as 

(63) B(z) = J B z (A)e^dA 

where B Z (X) is the Fourier component of a field perturbation with wavelength A. The 
total phase shift from magnetic field inhomogeneities is 

(64) A<j> B = J T Bz (X)B z (X)d\ 

Here Tg z (A) is the interferometer's magnetic inhomogeneity response function. 

As with gravity above, we compute the differential response function ATgz(A) = 
\{T B z)s7 - {T B z)$5\ between 87 Rb and 85 Rb (see Fig. 13). The differential response 
arises because the isotopes have different second order Zeeman coefficients a, as well as 
different masses. This sensitivity curve drives our magnetic shield design requirements, 
as discussed in Section 5 '3. 3. 

5'3. Controlling potential systematic errors. - 
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5'3.1. Rotation of the Earth. The largest systematic term in the phase shift expansion 
for a dual species differential interferometer after propagation reversal is due to the 
rotation of the Earth. Specifically, a differential acceleration due to the coriolis force 
occurs if the isotopes have different transverse velocities Sv x (Table II, term 2). Reducing 
this phase shift below our systematic threshold would require 5v x < 1CP 11 m/s, which is 
challenging. However, this specification can be relaxed by artificially making the rotation 
rate zero. 

To a good approximation^), the atoms are only affected by the Earth's rotation 
through their coupling to the laser, so the coriolis acceleration can be eliminated by 
rotating the laser in the opposite direction of the Earth's rotation. In order to calculate 
the effect of this rotation compensation, we performed the phase shift calculation using a 
rotating k c ff. Following the work of [14], we use a retro-reflection configuration to deliver 
the laser beams ki and k 2 to the atoms. We rotate k e ff by actuating the retro-reflection 
mirror. As a result, the incoming beam remains pointing along the z-direction and only 
the reflected beam rotates. With this configuration, k e ff is given by 

(65) k off = -2fcn(n • k) 

where n is the time-dependent unit normal vector of the retro-reflection mirror, and k 
is a unit vector in the direction of the fixed delivery beam. Notice that the direction of 
k ff rotates as desired, but its length now depends on anglc( 4 ). 

The resulting phase shift list appears in Table III, with SQ y and S£l z the errors in the 
applied counter-rotation rate. Assuming a transverse velocity difference of Sv x ~ 1 /zm/s, 
these rotation compensation errors must be kept below 10~ 5 riEarth ~ 1 nrad/s. Methods 
for measuring angles with nanoradian precision have already been demonstrated [36] . In 
order to actuate the mirror at this level of precision we can use commercially available 
sub-nm accurate piezo-electric actuators along with active feedback. 

Notice that not all rotation-related phase errors are removed by rotation compensa- 
tion. Terms that arise from the differential centrifugal acceleration between the isotopes 
(e.g., Table III terms 3 and 4) are not suppressed. Physically, this is a consequence of the 
fact that the retro-reflection mirror that we use to change the laser's angle is displaced 
from the center of rotation of the Earth by R e . Therefore, although we can compen- 
sate for the angle of the laser by counter-rotating, the retro-reflection mirror remains 
attached to the rotating Earth, leading to a centrifugal acceleration of the phase fronts. 



( 3 ) The atoms are also weakly coupled electromagnetically and gravitationally to the local 
environment, which is fixed to the rotating Earth. These cross-couplings to rotation are generally 
not important because the dominant gravitational interaction with the Earth is spherically 
symmetric, and all electromagnetic interactions with the atom (e.g. with the applied magnetic 
bias field) are naturally small. 

( 4 ) This small change in magnitude of k e a does not lead to any problematic phase errors in the 
interferometer since the total angle through which the laser rotates is only A9 = 20Earth7 1 ~ 
1(T 4 rad, and the effect is 0(A6> 2 ). 
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Phase shift Size (rad) Fractional size 



2 -k ctf T zz Sv z T 3 

3 -U- ~) hk 2 cS {l 2 v T 3 

4 -~3k eS Q. 2 y 6v z T 3 

5 -k cB T zz SzT 2 

8 -(^-^)^ cff B (d z B)T 2 

10 -2fc eff< 5 t ; ;c 50j / r 2 



1.66 x IO -2 


5.83 


x 


lO" 11 


1.44 x 10" 4 


5.05 


x 


io- 12 


-5.40 x IO -5 


1.90 


x 


10 -13 


-4.68 x IO -6 


1.64 


x 


IO" 14 


8.93 x IO -7 


3.14 


x 


IO" 15 


-1.61 x IO -7 


5.65 


x 


10 -16 


1.03 x 10" 7 


3.63 


x 


10 -16 


-9.94 x IO -8 


3.49 


x 


10 -16 


3.46 x IO -8 


1.22 


x 


10 -16 


-3.35 x IO" 8 


1.18 


x 


10 -16 



Table III. - Differential phase shift list with rotation compensation. Terms 1, 3, 6, and 7 will 
be suppressed by the propagation reversal technique described in Section 5'2. 



After propagation reversal, the only term of this type that is significant is ~ k c sQySv z T 3 
(Table III term 4). However, this term is smaller than and has the same scaling with 
experimental control parameters as the gravity gradient phase shift (Table III term 2), 
so the constraints described in Section 5'3.2 to suppress the gravity gradient terms are 
sufficient to control this centrifugal term as well. 

One potential obstacle in achieving the required transverse velocity constraint of 5v x ~ 
1 /im/s is the expected micro-motion the atoms experience in the TOP magnetic trap 
prior to launch [37] . Micro- motion orbital velocities in a tight TOP trap such as ours can 
approach <~ 1 cm/s in the transverse plane. Although the differential orbital velocities 
are suppressed by the 87 Rb- 85 Rb mass ratio, the resulting Sv x ~ 100 ^m/s is still too 
large. This problem can potentially be solved by adiabatically reducing the magnetic 
field gradient and increasing the rotating field frequency prior to launch. 

5'3.2. Gravity gradients. The largest systematic background after rotation com- 
pensation is due to the gravity gradient along the vertical direction of the apparatus 
(T zz — d z g z ). Since gravity is not uniform, the two isotopes experience a different av- 
erage acceleration if their trajectories are not identical. This effect causes a differential 
phase shift proportional to the initial spatial separation and initial velocity difference 
between the isotopes (see Table III, terms 2 and 5). Assuming a spherical Earth model, 
the gravity gradient felt by the atoms is T zz ~ 3 x 10~ 16 .g/nm, which means that the 
initial vertical position difference between the isotopes Sz must be < 1 nm and the initial 
vertical velocity difference Sv z must be < 1 nm/s in order to reduce the systematic phase 
shift beneath our threshold. 

The experiment is designed to initially co-locate the two isotope clouds at the nm 
level by evaporative cooling both species in the same magnetic trap. For trapping, we 
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use the state |F = 2, rap = 2) for 87 Rb and \F = 3, rap = 3) for 85 Rb since these states 
have the same magnetic moment [38] . The mass difference between the isotopes leads to 
a differential trap offset in the combined magnetic and gravitational potential given by 

(66) Az trap = — - 

where Am is the 87 Rb- 85 Rb mass difference, B" is the magnetic field curvature of the 
trap, and /is is the Bohr magneton. Our TOP magnetic trap is designed to provide a field 
curvature B" ~ 4 x 10 5 Gauss/cm 2 which reduces the trap offset to Az trap w 10 nm. The 
resulting systematic error is <~ 10~ 14 g, but it can be subtracted during the analysis given 
a knowledge of Aztrap at the ~ 10% level. This offset can be inferred from a measurement 
of the field curvature B" of the trap (e.g., by measuring the trap oscillation frequency). 
The gravity gradient must also be known, but this can be characterized in situ by using 
the interferometer as a gradiometer [15]. 

Control of the gravity gradient phase shift (Tabic 111 term 2) requires that the dif- 
ferential launch velocity be Sv z < 1 nm/s. Therefore we cannot employ standard launch 
techniques (e.g., moving molasses) since the velocity uncertainty is fundamentally limited 
by the photon recoil velocity vpi <~ 6 mm/s due to spontaneous emission. Instead, the 
atoms are launched using an accelerated optical lattice potential [27] . We launch the two 
species using the same far-detuned (~ 200GHz) optical lattice, coherently transferring 
<~ 2200hk of momentum to each cloud. Because the two species have different masses, 
they have different Bloch oscillation times tr = , where a is the lattice acceleration. 

J ^ ma 7 

As a result, after adiabatically ramping down the lattice potential, the two species are in 
different momentum eigenstates since they have absorbed a different number of photons. 
The differential velocity after launch is then 

(67) Sv L = hk cS f^^^l) 

\m S5 m s7 J 

where Ng 5 and N 87 are the number of photons transferred to 85 Rb and 87 Rb, respectively. 
We choose the integers Nss = 2168 and Ns7 = 2219 such that their ratio is as close to 
the isotope mass ratio as possible, resulting in Svp ~ 12 /im/s. After launch, we can 
perform a velocity selective transition to pick out a common class from the overlapping 
distributions of the two isotopes, which at the expense of atom number could conceivably 
allow us to achieve our differential velocity constraint. 

There are several additional ways to reduce the gravity gradient systematics beyond 
precise control of the launch kinematics. We can implement a 4-pulse sequence (-| — 7r — 
7r — J) which suppresses all phase shift terms oc T 3 at the cost of an order one loss in 
acceleration sensitivity [5] . This eliminates the velocity dependent gravity gradient phase 
shifts but would still require that we maintain tight control over the initial differential 
vertical position between the isotope clouds. Secondly, we can potentially reduce the 
local gravity gradient T zz by applying appropriate trim masses around the apparatus. 
It has been shown [34] that in principle a local mass distribution can effectively cancel 
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the gravity gradient of the Earth for a 10 m-scale apparatus. Reducing T zz by an order 
of magnitude would relax our initial position constraint to the level provided by the 
expected value of A2; trap , thereby removing the requirement for subtraction during data 
analysis. 

5'3.3. Magnetic fields. The magnetic field phase shift appearing in Table III (term 
8) constrains the maximum allowed linear field gradient to d z B < 0.1 nT/meter. In 
the interferometer region, the measured gradient of the Earth's field is ~ 3 /iT/m, and 
therefore we require a shielding ratio of at least <~ 5 x 10 4 . In addition to suppressing the 
field gradient, the magnetic response function (see Fig. 13) indicates that the field must 
be uniform on length scales ~ 1 cm. Large magnetic shields with similar performance 
have been demonstrated [39]. The magnetic shielding for our interferometer region is 
provided by a three-layer concentric cylindrical shield made of high permeability mate- 
rial. To maintain a pristine magnetic environment, we use an aluminum vacuum chamber 
and non-magnetic materials inside the shielded region. 

In order to verify the performance of the magnetic shield, we must characterize the 
field. As with gravity inhomogencities, the atom interferometer can be used to map the 
local magnetic field in situ, in this case by using a magnetic field sensitive (mp ^ 0) 
state [40]. 

6. Conclusion 

In these notes we have given an overview of the light-pulse method, and discussed 
applications to inertial navigation and a test of the Equivalence Principle. As the field 
continues to progress, we see two trends evolving. First, a steady evolution in the tech- 
nology associated with laser manipulation of atoms will lead to progressive development 
of compact, field-ready inertial sensors. Such developments include integrated photonics 
packages for the laser and optics paths, and integrated optics bench and vacuum systems. 
On the other hand, continued evolution of high performance science instruments will ex- 
tend the physics reach of this technology. For example, next generation light-pulse atom 
interference systems appear capable of placing superb limits on atom charge neutrality 
[41], making terrestrial tests of General Relativity [34, 42], and possibly detecting low 
frequency gravitational waves [43, 44]. 
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